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SUMMARY 

The flowfield characteristics in rocket engine coolant channels are analyzed by means of a numerical model. 
The channels are characterized by large length to diameter ratios, high Reynolds numbers, and asymmetrical heating. 
At representative flow conditions, the channel length is approximately twice the hydraulic entrance length so that fully 
developed conditions would be reached for a constant property fluid. For the supercritical hydrogen that is used as the 
coolant, the strong property variations create significant secondary flows in the cross-plane which have a major influ- 
ence on the flow and the resulting heat transfer. Comparison of constant and variable property solutions show substan- 
tial differences. In addition, the property variations prevent fully developed flow. The density variation accelerates the 
fluid in the channels increasing the pressure drop without an accompanying increase in heat flux. Analyses of the inlet 
configuration suggest that side entry from a manifold can affect the development of the velocity profile because of vor- 
tices generated as the flow enters the channel.Current work is focused on studying the effects of channel bifurcation 
on the flow field and the heat transfer characteristics. 

TECHNICAL DISCUSSION 

Regenerative cooling is normally used in rocket engine combustors to protect the walls from the severe ther- 
mal environment of the combustion process. In most cases, the coolant passages are rectangular in shape and run the 
length of the combustor with the coolant entering at the supersonic end, passing through the throat region, and exiting 
near the injector face. In cryogenic engines, the coolant is generally liquid hydrogen, and the elevated pressures are 
sufficient to ensure that the fluid remains supercritical over the entire length of the cooling passage. 

Although the design of these coolant passages is critical to the operation and life of the engine, their heat 
transfer characteristics are only poorly understood. To maximize our understanding of the fluid dynamics and the as- 
sociated heat transfer characteristics of the coolant passages, the limiting case of a straight rectangular channel (which 
is currently being tested at NASA LeRC) is studied. Even with straight geometries, the flowfield inside the tubes re- 
mains quite complex because of the asymmetrical heating, the fin effect, the presence of a supercritical coolant with 
its characteristically strong property variations and the high Reynolds numbers. Because of the complex three-dimen- 
sional flow patterns in these coolant passages, it is anticipated that an optimum channel aspect ratio exists, and a vali- 
dated CFD code is an important tool in assessing where this optimum occurs under different conditions. 

At supercritical conditions all fluid properties (density, viscosity, thermal conductivity and specific heat) be- 
come strong functions of both pressure and temperature [1]. These variations are particularly strong near the critical 
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point. Temperature dependence is especially significant in the present problem because of the strong temperature gra- 
dients near the hot wall of the coolant passage. 

In the presence of property variations, the Navier-Stokes equations and the energy equation are coupled and 
must be solved simultaneously. The corresponding system of governing equations is: 

r s + | (E - Ev)+ li <F - Fv)+ ^ (G - Gv) =° 

Q=(P,u,v,w,T) T is the vector of dependent variables;(£, t|, Q represents a general non-orthogonal coordinate system; 
Tis a pre-conditioning matrix; E, F, and G are the conservative flux vectors while Ey, F v and G v are the viscous flux 
vectors [2], The equation set is completed by a tabular equation of state relating the density to the local temperature 
and pressure. Additionally, auxiliary equations for the pressure and temperature dependence of the viscosity, thermal 
conductivity and specific heat are also specified. The Baldwin-Lomax algebraic turbulence model [3] is used to accom- 
modate turbulent flow. The turbulent Prandtl number is chosen as 0.9 throughout the calculation. 

Although our interest is in steady flow, the equations of motion are written in unsteady form for computational 
purposes and are solved by a four stage Runge-Kutta explicit scheme [4]. All derivatives of the Navier-Stokes equa- 
tions are represented by centered differences. Local time stepping is used to achieve maximum convergence and a 
small amount of fourth order artificial dissipation is added to prevent odd-even splitting. Since the time derivatives dis- 
appear in the steady state, some license may be taken with them to counter the stiffness caused by the low Mach num- 
bers in the supercritical fluid. This stiffness is circumvented using an artificial compressibility expression [5,6]. During 
the process of convergence, all the time derivatives approach zero and the numerical solution satisfies the proper steady 
state Navier-Stokes equations. 

As boundary conditions for the present calculations, we specify the velocity and temperature profiles at the 
inlet, and uniform pressure at the exit plane. The remaining flow variables at the inlet and outlet are computed from 
the Riemann variables determined from the method of characteristics [4]. The velocity boundary conditions on the pas- 
sage wall are the traditional no-slip conditions. These are augmented by enforcing the normal momentum equation to 
obtain the wall pressure. 

Heat transfer conditions representative of those in combustor coolant passages are simulated by selecting ap- 
propriate Dirichlet (temperature) boundary conditions around the periphery of the duct. Uniform temperature condi- 
tions are assigned on both the inner (combustor) and outer (ambient) sides of the passage while a linear distribution of 
temperature is used on the side walls to simulate the fin effect. The inlet temperature is specified as 40 K. Future plans 
are to couple the present fluids solution with a heat conduction formulation for the combustor walls to obtain the cou- 
pled solution. 

The straight duct study was conducted to give an understanding of the fundamental physical effects of variable 
properties, turbulence modeling, high length to diameter ratios and inlet conditions. The next step in the study is to 
investigate the effects of channel bifurcation on the flow and temperature fields. A representative geometry is shown 
in Figure 1 . 

RESULTS 

The geometry and flow conditions for this study were chosen to parallel an experimental study currently being 
conducted at NASA LeRC [7], The length of the experimental channel is 137 mm which corresponds to 1 12 hydraulic 
diameters. The Reynolds number based on the inlet conditions and hydraulic diameter is 500,000. For the current study 
we also studied channels of shorter lengths and flow at the lower Reynolds number of 10,000. Computations have also 
been completed for both constant and variable (supercritical hydrogen) property flows. 
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One of the important issues in the straight tube studies is the effect of variable properties on the pressure drop 
and heat transfer characteristics of the channels. The property variation of super-critical hydrogen significandy alters 
the flow field characteristics and it is important to assess its effect on the channel cooling. To evaluate the nature of 
these differences a comparison of the flowfields based upon constant and variable properties was conducted. This com- 
parison was made at a Reynolds number of 10,000. 

The mid-plane, streamwise direction contour plots of the pressure, the axial velocity and the temperature vari- 
ations are shown in Figure 2. The constant property results are shown on the right while the variable property results 
are shown on the left. The cross-stream direction has been stretched by a factor of five in theses figures. As can be seen, 
the variable properties cause the velocity profiles to be very non-uniform across the tube, and in addition, lead to an 
acceleration by a factor of about two. The higher flow speeds cause a lesser, but noticeable effect on the temperature. 

A comparison of the flow fields for the two different Reynolds numbers 10,000 and 500,000 shows the viscous 
and thermal boundary layers were significantly thinner at the higher Reynolds number. Despite the fact that the thermal 
effect are limited to a small portion of the flowfield there is still an acceleration of the flow due to the density variations 
although not as significantly as in the lower Reynolds number. 

Another issue studied was the effect of inlet velocity profiles on the downstream fluid dynamics and heat 
transfer. The inlet conditions studied were slug flow, fully developed velocity profiles and entry through a 90 degree 
bend (constant property). The first two represents a bounding study between two extremes in inlet profiles. The last 
condition represents an approximation to the type of entry that is used in the experimental configuration. A comparison 
of fully developed velocity profile with the slug flow velocity profile shows that the inlet condition has a measurable 
effect on the velocity contours in the exit plane (Fig. 3). The contours in the exit plane are qualitatively similar, but the 
results for the slug flow inlet show a peak velocity of 4.0 m/s in the exit plane whereas those for the fully developed 
inlet show a peak of 3.3 m/s. The corresponding pressure drop for the slug flow inlet is about 15 percent larger than 
that for the fully developed inlet profile. By contrast, the fluid temperatures at the exit plane for both cases are about 
the same. Both calculations show that the fluid has been heated to the point where the minimum temperature in the exit 
plane is 45 K, although more of the fluid appears to have been heated in the slug flow inlet case. In agreement with this 
latter observation, computation of the heat flux through the walls shows that the heat addition is larger (by about 15 
percent) for the slug flow inlet case. These differences between the two calculations are not major, but they do suggest 
that care should be taken in defining the inlet conditions to be sure to obtain maximum accuracy in the coolant channel 
predictions. 

The results for the second inlet condition, the inlet with a 90 degree bend are presented in Figure 4. The ve- 
locity contours are quite similar in magnitude, although the profiles for the 90 degree bend inlet are slightly asymmet- 
ric. Similar comparison of the temperature profiles shows the presence of the bend in the inlet increases the minimum 
temperature from 90 K to 1 10 K, suggesting that the vortex created by the turn at the inlet has some effect on the heat 
flux. 

Comparisons of the streamwise development of the various flow characteristics are given on Figs. 5 and 1. 
The pressure drop predictions for the two cases is nearly identical (not counting the pressure drop across the inlet bend 
itself)* suggesting that the strength of the vortices developed by the bend in the inlet was not sufficient to affect the 
downstream pressure gradient. The velocity profile in the bent inlet calculation shows a dramatic thickening of the 
boundary layer downstream of the inside portion of the curved inlet. This arises from a flow separation near the bend. 
The presence of flow separation and the vortex generated by the curved inlet suggests that the type of connection be- 
tween the manifold and the coolant channels can have an effect on the downstream heat transfer. 
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Figure 1: Representative sketch of channel bifurcation geometry. 
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Figure 2: Steamwise contours of the pressure, velocity, and temperature on the centerplane for an L/D h of 28. Vari- 
able Properties: P^ = 20., P max = 220.; = 0.5, u max = 4.0; T ^ = 70, T max = 600 (top wall). Constant Properties: 

^inc = F max — 42 u inc = 0.17, u max =1.4, T mc = 70, T max = 600 (top wall). 
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Figure 3: Contour of axial velocity and temperature at 

the exit plane (L/D h = 28) for slug flow and fully devel- 
oped inlet with variable properties. Slug: u^ = 0.4, u,^ 
= 4.0, Tine = 70, Tjnin = 45. Fully Developed: u^ = 0.4, 
u max = 3.3; Tj Jlc = 70, T max = 600 (right wall). 


Slug Flow Inlet ,, Bend Inlet 

Velocity 


Temperature 



Figure 4: Contour of axial velocity and temperature in 

the exit plane (L/E> h =28) for slug flow and bend entrance 
inlet with constant properties. Slug: u^ = 0.4, u max = 4.0, 
T inc = 70, Tmin = 45. Bend Inlet: u mc = 0.4, u max = 3.3; 
Tjnc = 70, T max = 600 (right wall). 





Figure 5: Effects of bend on flowfield (constant properties). There are 28 hydraulic diameters after the bend. 

Pine = 0.5, P max = 5; u mc = 0.1, umin = -0.1, u max = 1 .5; T mc = 70, T max = 600 (at bottom wall), Tmin = 40 (at inlet). 
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